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Abstract 

Aiming at a microscopic description of heavy ion collisions in the beam energy region of about 
10 A GeV, we extend the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model by in- 
cluding a relativistic mean field, in-medium baryon-baryon cross sections and three-body collisions. 
The model is then compared with experimental data for central Au+Au collisions at 2-10 A GeV 
and central Fb+Fb collisions at 30 and 40 A GeV on the proton rapidity spectra, the midrapidity 
yields of 7r + , K and (A + E°), and the transverse mass spectra of 7r ± and K^. The three-body 
collisions increase the inverse slope parameters of the hadron mj_-spectra to a good agreement with 
the data. 

PACS numbers: 24.10.Lx; 24.10.Jv; 25.75.-q; 25.75.Dw 
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I. INTRODUCTION 



The Boltzmann-Uehling-Uhlenbeck (BUU) transport theory, based on binary collisions 
and on the propagation of particles in a selfconsistent mean field, is a very useful tool in 
understanding the heavy ion collisions (HIC) in the energy region from the Fermi energies 
(Eu h ~ 30 A MeV) up to the relativistic energies (E uh ~ 2 A GeV) (c.f. Refs. [l 3 3 L 4^5, 
y, |7| for the description of the model in various numerical realizations and Refs. [8J, |9|, lid. Ill] 
for the results below 2 A GeV). 

There are, however, systematic deviations from experimental data on pion and kaon 
production at E\ ah > 2 A GeV [l2, Q, Q: the pion multiplicity is systematically overpre- 
dicted, while the slopes of the K + transverse mass spectra are too steep in the transport 
calculations. It has been advocated in [l3j, that the too soft kaon m_i_-spectra are caused by 
missing the formation of a nonhadronic phase, which should create an additional pressure 
accelerating the kaons. On the other hand, the three-fluid hydrodynamical calculations [ijj] 
have apparently been quite successfull in reproducing the m^-spectra of hadrons at AGS to 
SPS energies by using a hadronic equation of state. This points to the idea of not enough 
thermalization produced by the microscopic transport models rather than to the importance 
of the nonhadronic degrees of freedom. 

The many-body collisions — usually missed in the current transport calculations at high 
energies — could serve as an additional source of the thermalization. The role of many-body 
collisions grows with baryon density, which reaches values of order of 1-2 fm~ 3 at the energy 
region of 10-20 A GeV [161 ] . These are typical energies of the future Compressed Baryonic 
Matter experiment at the Facility for Antiproton and Ion Research in Darmstadt ~v\ . 



A simple estimate of the gas parameter 



181 ] at the baryon density p# = 10po, which is 



the maximum density reached in a central Au+Au collision at 20 A GeV with p = 0.17 
fm -3 being the normal nuclear matter density, is 



7gas 



{a/nf' 2 p B ~ 2 , (1) 



where a ~ 40 mb is the asymptotic high-energy value of the total pp cross section in vacuum. 
In thermally equilibrated nuclear matter, one can neglect the Lorentz contraction of the 
interaction volume |7r(cr/7r) 3 / 2 (c.f. Eqs. flA2l) . fT4T)j) and Fig. [2] below). Since 7 gas > 1, one 
concludes, that the applicability condition of the binary collision approximation is violated 
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fcf. Ref. 



19l|). As we have recently shown [20J, at the maximum compression stage of the 



central Au+Au collision at 20 A GeV, iV-body collisions with N > 6 should dominate. 

One way to describe this complex physical situation could be, indeed, an introduction of 
the new degrees of freedom, e.g. like the formation of a quark-gluon plasma. This is quite 
a challenging problem for the microscopic transport theories, which is, however, beyond the 
scope of our present study. Another way is to simulate the many-body collisions within 
a theory containing only the hadronic degrees of freedom. At beam energies below 1-2 



A GeV, transport models taking into account A-body col 



constructed by several authors (c.f. Refs. 



21 



22 



23 



24 



isions (with N > 3) have been 
251]) . The difficulty, which appears 



at high baryon densities reached in the 10 A GeV region is that the gas parameter is not 
small, and, therefore, such a theory can not be formulated with the vacuum cross sections 
as a decomposition of the collision integral in a series of powers of 7 gas (c.f. Ref. |18|). 
In other words, the vacuum cross sections should be screened at high densities by particles 
surrounding the colliding pair 26j. The screening effect appears also as a consequence of the 



Dirac-Brueckner calculations (c.f. Refs. 



27 



28] ), where the nuclear mean field and the in- 



medium reduced cross sections are derived from the same fundamental vacuum interaction. 
Using in-medium reduced cross sections would reduce the relative contribution of the many- 
body collisions, since the gas parameter will be smaller. This would give a more solid ground 
,oJ;he kinetic theory, which is based on the small parameter 7 gas . The Fermi liquid theory 



29, 



30] gives a similar picture. There, the liquid of the real particles is equivalently replaced 



by the gas of quasiparticles. 

In the present work, we develop a transport model which contains a Walecka-type baryonic 
mean field and the in-medium reduced baryon-baryon cross sections. The model is then 
extended by including three-body collisions. The mean field and the in-medium reduced 
cross sections lead to less thermalization, while the three-body collisions counterbalance 
this effect. We will show that the three-body collisions raise the inverse slope parameters of 
the hadron spectra in the central HIC at beam energies of 2-40 A GeV. In particular, the 
measured K + transverse mass spectra are well described by calculations with the three-body 
collisions. 

The structure of the paper is as follows. In Sect. II we briefly describe the GiBUU model 
j?! concentrating on its new ingredients: the relativistic mean field, the in-medium baryon- 
baryon cross sections and the three-body collisions. Sect. Ill contains numerical results. In 
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Sect. IV we summarize and discuss our results. 



II. GIBUU MODEL 



1™ 



Our calculations are based on the GiBUU model in a new version of Ref. [7|, written 
using FORTRAN 2003 in an object-oriented way. The model describes a nucleus-nucleus 
collision explicitly in time as a sequence of the baryon-baryon, meson-baryon and meson- 
meson collisions, baryon and meson resonance excitations and decays. Between the two- 
body collisions the particles propagate in a selfconsistent mean field. The baryon-baryon 
collisions at < 2.6 GeV are treated within the resonance model, while at a/s > 2.6 GeV 
the FRITIOF model is applied. For the meson-baryon collisions, the FRITIOF model is 
used at ^Js > 2.2 GeV. We applied the energy- dependent strangeness suppression factor 

0.4 for ^ < 5 GeV 

P(s) i 

7s = = { 0.433 - ^v/i [GeVH 1 for 5 GeV < y/s < 20 GeV (2) 

from Ref. 



150 

0.3 for v/i > 20 GeV 



3jj instead of the default FRITIOF value of 7 S = 0.3. In the meson-baryon 
collision case, the FRITIOF mechanism of the double-string excitation and decay has been 

n 

improved by adding the qq annihilation channel as described in Ref. [141 ] . 

The particles produced in string decays are not allowed to interact with their default 
cross sections up to some proper time interval, called formation time. We use a value of the 
formation time 77 = 0.8 fm/c for all baryons and mesons. If one of the colliding particles is 
still in the formation interval, we call it a prehadron. The total cross section of the prehadron 
interaction with another particles is scaled according to the constituent quark model (c.f. 



Ref. [32j]): If the prehadron interacts with a hadron, the total cross section is multiplied 
by the factor of iVi ea ding/3 (iV"i ea ding/2) for the pre-baryon (pre-meson), where iV"i ea din g is the 
number of leading (anti)quarks in the prehadron, i.e. the number of (anti)quarks which were 
existing in the parent colliding particles for the given prehadron. If both colliding particles 
are prehadrons, the total cross section of their interaction is multiplied by the product of 
the two corresponding factors. 

The rescaling factor of the prehadron cross sections is quite important for the description 
of the produced particle multiplicities. In previous GiBUU calculations of HIC at 2-40 A 
GeV [141 ]. this factor had been chosen to be 1/3 for the prehadrons containing at least one 



leading quark and zero otherwise. This explains the higher particle abundancies produced 
in the cascade mode in the present paper as compared to the results of Ref . 14| . Another 



distinction to the calculations of Ref. 



141 ] is that in the present work we miss the low energy 



(y/s < 2.6 GeV) baryon-baryon channel of kaon production BB — > BYK. However, this 
channel was found to be negligible for the beam energies above 4 A GeV 14]. 

In this work, we improve the previous GiBUU calculations [ljj] by implementing a rela- 
tivistic mean field (RMF), in- medium cross sections and three-body collisions in the model. 
The in-medium cross sections for the high energy baryon-baryon collisions computed ac- 
cording to the RMF model have already been introduced in 14j. However, the RMF has 
not been used for the particle propagation in |l4j. Below, we describe in-detail these new 
ingredients of our model. 



Relativistic mean field 



In distinction to the earlier GiBUU calculations [5|, |6|, |8|, |9|, |H|, [32 ] , where the nonrelativis- 
tic momentum- dependent potential of Ref. [sjj was used, in the present work we employ 
the relativistic treatment of the baryonic mean field. In spite of the Lorentz invariant imple- 
mentation of the potential from 33j , when it is calculated in the local rest frame of nuclear 
matter, the RMF is better suited for description of the high baryon densities reached in 
a central heavy ion collision at about 10 A GeV beam energy. Here, a consistent nuclear 
equation of state at high densities is of the primary importance. RMF models provide a 



quite good descrip 
neutron stars 



34 



ion of both the collective flow in HIC and of the high density interior of 



351 ]. They are also successfully applied in nuclear structure calculations 



361 ] . The one weakness of the RMF models, the too repulsive nucleon-nucleus interaction 



at high momenta (c.f. [34| and refs. therein) , is of minor relevance for the present results 
for particle production from high-density equilibrated nuclear matter. 



We use the relativistic mean field Lagrangian density in the form given in Ref. 



C = ijj[y(id - g^u) 



+ ^) 2 



U{a) 



(3) 



where ip is the nucleon field, a and uo are the scalar-isoscalar and the vector-isoscalar meson 
fields, respectively; F^ u = d^uj u — d u uj^. We neglect the isovector meson and the electromag- 



5 



netic contributions. The Lagrangian density (J3]) contains the nonlinear self-interactions of 
the cr-field: 

U(v) = \<° 2 + \^ + \g,° A ■ (4) 
The Lagrange's equations of motion for the nucleon, a- and u>-fields are written as follows: 

[j(id - g u u) - m nuc - g^a}^ = , (5) 
d^o + = -gM , (6) 

8^ + ml^ = g^ 1 ^ . (7) 

The a- and u;-fields are treated as classical ones, i.e. we replace these fields by their expec- 
tation values in Eqs.(j51 El E])- Assuming that the meson fields are varying much more slowly 
in time and space with respect to the nucleon field, we consider the plane-wave solutions of 
Eq.©: 

oc exp(=Fzpx) , (8) 

where x = (t, r), p = (p°, p), and the upper (lower) sign corresponds to the nucleon (antin- 
ucleon). The dispersion relation is then obtained from (jSJ): 



p° = ±g^ + yj(p T g u u ? + (m nuc ) 2 , (9) 

where 

m nuc = m nuc + 9o<T (10) 

is the nucleon effective (Dirac) mass. 

The distribution function f(x, p) of a given particle species in the phase space (r, p) is 
defined now such that f(x, p) 3 ^!^ = ( num ber of particles of that species in the phase 
space element d 3 rd 3 p), where g — 4 is the spin-isospin degeneracy. The space-time evolution 
of the (anti)nucleon phase space distribution function is described by the BUU equation 

df dp® df dp® df 

dt ~^ dp dr dr dp co11 ^] ' ( ) 



where — in spirit of the Fermi liquid theory [29|, l30( — the single-particle energy (jHJ) is used 
as a one-body Hamiltonian function. The l.h.s. of Eq. ffTT]) describes the propagation of 
particles in the mean field, while the r.h.s. is the collision integral. The explicit form of J co n 
for elastic scattering is given below (see Eqs. f|T5l) . fr22l) . flA4j) ). 
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It is convenient to perform a variable transformation in Eq. ljTT!) by using the kinetic 
four- momentum 

p^=p^T9^. (12) 

The distribution function f*(x, p*) in the kinetic phase space (r, p*) is defined such that 
f*(x, P*) gd ( 2 ^)3 P — (number of particles in the kinetic phase space element d 3 rd 3 p*). Since 
d 3 p = d 3 p*, which is valid for the momentum-independent to- field, one gets: 

f*(x,p*) = f(x,p) . (13) 

Expressing the l.h.s. of Eq.t jTTi) in terms of the kinetic quantities we obtain the transport 



equation (c.f. Ref. [341] ) : 



(Po) _1 K^ + (±9.P:F- + mi nc {d«<J)d P a ir(x, P*) = /ooutT] , (14) 



rmc/ 



where a = 1, 2, 3, and p$ is determined from the mass shell condition (p*) 2 = (m* 

The collision integral, generally, can be expanded in the number of colliding particles (c.f. 
Ref. [22j and refs. therein): 

-fcoll[/*] = icoll,2b[/*] + -^coll,3b[/*] + -^coll,4b[/*] + ••• ■ (15) 



We will restrict ourselves to the first two terms in ([15]) only. It is straightforward to write 
down the elastic contributions to J co ii,2b an d I C oii,3b m the case of identical fermions: 
reiasticr^i 1 ™i / gd 3 p\,m\ l f gd 3 p\ l m1 1 , f gd 3 p\m\ 



relastic r r*~ 
i coll,2bL/l. 



2! pf J (2n) 3 pf J (2n) 3 p*? J (2vr) 3 ^ 



x|M|? 2 _ 1/2 , (2n)U^(p 1 +p 2 -p v - p v ) - ft ft ft %) , (16) 

relastic r f*' 
i coll,3bUl. 



L;sil 1 m\ f gd 3 p\ l m\ l f gd 3 p\,m\, f gd 3 p* 3/ m* 3 , 



3\2\ pf J (2vr) 3 p*P J (2vr) 3 p*? J (2vr) 3 _p*? 

/ Mil / W ^ + »» + » - »* - »» - *-) 

xUifififtnit - nnnnJifi) , 



x 



where |M|f 2 _ > i'2' an d |^|i23-»i'2'3' are the spin- and isospin-averaged over initial and fi- 
nal states matrix elements squared for two- and three-body collisions, respectively; /* = 
f*(x, pf), f* = 1 — /* (i = 1, 2, 3, 1', 2', 3'). The matrix elements satisfy the detailed balance 
relations 



= |M|^ 1/2 , , (18) 



l^1l'2'3'-»123 ~~ I^1l23->1'2'3' ' (19) 



which are used in Eqs. fll6l) . (|T7|) . The normalization of Bjorken and Drell [37j is used for 
the matrix elements, which leads to the appearance of the fermion Dirac masses in Eqs. 
(IT6|) . (jl7p . In the case of the momentum-independent a-field used in this work, we have for 
the identical fermions m* = m 2 = = m\, = m 2 , = m^,. However, in Eqs. (jl6p . (|17|) and 
everywhere below in the expressions for the collision integrals, we keep, for clarity, different 
subscripts for the Dirac masses of different particles. 

Introducing the in-medium elastic differential scattering cross section 



da\ 2 ^ V2 , = (27r) 4 6 (4 \p 1 + P2 -p 1 ,-p 2/ )\M\l 
m\m 2 gd 3 p\,m\, gd 3 p 2 ,m 2 , 1 



ijfe (2tt) 3 ^? (2tt) 3 p*° 2! 



with 



(20) 



i* = J{p\p* 2 Y - (mimt) 2 (21) 



being the in-medium flux factor, one can rewrite the two-body elastic collision term as 
follows: 

gd 3 p* 2 



Cf 2 M =Jj^J d °^V2> v*M>tiftft ~ Ktiftfi,) , (22) 



where 



iy(pf P f) (23) 



U 12 

is the relative velocity of colliding particles. 

In agreement with the low-energy in-medium calculations of the elastic NN scattering 
cross section (c.f. Ref. [381), one can neglect the medium dependence of the matrix element 



|M|f 2 _ > i'2' approximately. This approximation will be used in the next subsection in order 
to evaluate the in-medium baryon-baryon cross sections (c.f. Eq. (1351) below). 

To solve Eq. (IHl) . the distribution function in the kinetic phase space is projected onto 
test particles: 

f27r) 3 AN 

n*, P*)~E *(r - r,(t))«y(p* - pt(t)) , (24) 

where A is the number of nucleons and iV is the number of test particles per nucleon. The 
centroids of the 5-functions are evolving in time between the two- or three-body collisions 
according to the following equations, which can be obtained by substituting (12111 into (JT41 



and putting the collision term equal to zero (c.f. Refs. [34j, |39(): 

dvi 



dt p 



(25) 



i Piu zpau | m nuc 30^* 

aL Pi Pi 

with a = 1, 2, 3 and z/ = 0, 1, 2, 3. 

The a- and cu-fields are calculated from the Klein-Gordon equations (jSJ),© by neglecting 
the derivatives of the fields in space and time: 

m 2 a a + g 2 cr 2 + g 3 v 3 = -gaps , (27) 
mlu v = g u f B . (28) 

(29) 

The scalar density p s (x) =< ip(x)ip(x) > and the baryon current j B (x) =< ip(x)^ p ip(x) > 
are expressed via the (anti)nucleon phase space distribution functions: 

Ps( X ) — J {2tc) 3 P*^ C ^ nuclcon ^* E ' P ) + /antinucloon( :r ; P )) ) (30) 
Jb( X ) = J "^^p^0"(/nucleon( a '5 P ) — /antinucleon ( X ! P )) ■ (31) 

Since we use the distribution functions in the kinetic phase space, Eqs.f l27|) . fl30l do not de- 
pend explicitly on the vector field to. This simplifies the selfconsistent numerical calculation 
of the meson mean fields strongly. 

We assume, for simplicity, the same coupling constants of all other baryons with the 
a- and a;-fields as for the nucleon. Correspondingly, in actual calculations of the scalar 
density and of the baryon current, the partial contributions from all the baryons present 
in the system are taken into account. The baryon test particles are propagated according 
to Eqs. fl25|) . ff26l) with a replacement m* uc — > m B = m B + g a a, where mg is the vacuum 
mass of the baryon B. The potentials acting on mesons are neglected. Thus, the mesons 
are propagated freely between the two- or three-body collisions. 

The numerical values of the RMF parameters were chosen according to the parameter set 
NL2 from Ref. [39[: m a = 550.5 MeV, m u = 783.3 MeV, g a = 8.50, g u = 7.54, g 2 = -50.37 



fm x , #3 = —6.26. This parameter set 40j produces the incompressibility K = 210 MeV and 



the ratio m* uc /m nuc = 0.83 at normal nuclear matter density. At high densities, the NL2 



34J]. This agrees with the BEVALAC 



parameter set gives a rather soft equation of state 
data on collective flow in HIC at about 1 A GeV beam energy 34]. 

The ex-field and the baryon four-current have been computed on the space grid with 
the cell sizes ~ (1 x 1 x 1/7) fm covering the collision zone. The smaller cell size in the 



longitudinal direction is needed in order to resolve the density profiles of the colliding nuclei, 
which are Lorentz contracted by the 7-factor in the center-of-mass (cm.) frame. We have 

n 

used the parallel ensemble technique (c.f. Ref. with N = 200 test particles per nucleon. 
To get smooth meson mean fields, the 5-function in coordinate space 5(r — i"i(t)) in Eq. fliH|) 
has been replaced by the Lorentz contracted gaussian 

Pl[r} ~ (2ir)*W ex Pj ^2 2# J (d j 

with L ~ 1 fm. The equations of motion fl25l) . fl26l) for the test particles have been solved by 
using the 0(At 2 ) predictor-corrector method. The space derivatives in the r.h.s. of Eq. (!26|) 
have been computed using the central differences, which produces a second order accuracy 
also in space. The numerical scheme conserves energy with an accuracy better than 3% of 
the cm. kinetic energy for the studied reactions. 

B. In- medium cross sections 

The most important part of the GiBUU model is the collision integral which includes 
the hadron-hadron cross sections. In the case of cascade calculation, neglecting any mean 
field effects, the vacuum cross sections are invoked. These cross sections are based either 
on the resonance model or on the phenomenological parameterizations of the experimental 
data with some reasonable extrapolations to the not-measurable channels (e.g. in the case of 
meson-meson collisions). The detailed description of the low energy resonance cross sections 

n n n 

is given m Refs. [3|,15|,16|. The high energy cross section parameterizations — used for the 



FRITIOF event generator — are explained in Ref. 32J. For the meson-meson cross sections, 
we refer the reader to Ref. |l4]. 

Using the RMF model in particle propagation requires also to introduce in-medium mod- 
ifications of the cross sections. This should be already clear, since the particle production 
thresholds include now the Dirac masses instead of the vacuum masses. To evaluate the 
cross sections in the case of calculations with RMF, we will apply two different schemes. 

In the first scheme, a so-called corrected invariant energy of the two colliding particles 1 
and 2 is computed as 

\/ S corr = Vs* - (m\ - 7TT,l) - (m% - m 2 ) , (33) 
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where 

S* = (pl + Pt) 2 ■ (34) 

The corrected invariant energy is an analog of the vacuum invariant energy, since the scalar 
selfenergies of the colliding particles are subtracted in the r.h.s. of Eq.(l33l. The quantity 
^/Scorr is then used in calculation of any reaction cross section <7i2-^x(\/s corr ) ■ Due to the same 
scalar selfenergies for all the baryons, this scheme ensures the correct threshold conditions 
for all binary processes except for BB production/annihilation. The last processes are, 
however, not important at the beam energies considered in the present work. 

In the second scheme, we follow the approach of Ref. [ijj]. There the in-medium cross 
section of a process B1B2 — > B^B^M^M^...M^ with B\,i and B3 4 as incoming and outgoing 
baryons and M 5j iA r as produced mesons, is expressed in a form: 

a med (V^) = Fa vac (^;) . (35) 

The modification factor F is 

_ m\m\m%ml I $^-2(^/3*; m%, m%, m* N ) 
mim 2 m z m A I* $at- 2 (a/s^; m 3 , m 4 , ...,m N ) ' 

where 

$ n (M;m 1 ,m 2 ,...,m n ) = J d^ n (P;p 1 ,p 2 , -,p n ) (37) 
is the n-body phase space volume with mf = p 2 (i = 1, 2, n) and M 2 = P 2 , 

d^(P; Pl ,p 2 , ...,p n ) = 8^{P-p x -p 2 - ... - Pn) (27r)32p o (27r)32p o • • • 7^)^ (38) 



is the element of the n-body phase space volume (c.f. Ref. [41]) : 



In = ?(V ' Scott, m 1 ,m 2 )y/ '-Wr , (39) 

J* = 9 (VF,mJ,m^)Vs* (40) 



are the vacuum and in-medium flux factors (c.f. Eq. (|2ip ) with 



q(^/s, mi,m 2 ) = \J (s + m 2 — m 2 ) 2 / (As) — m 2 (41) 

being the center-of-mass (cm.) momentum. Eqs.( l35"l) . (l3"6l) take into account the in-medium 
modification of the Dirac plane wave normalization (a factor of m\m\m\m\l (yn\m 2 mzm^\ 
of the flux factor and of the phase space volume. However, it is assumed that the matrix 

11 



element of the reaction is not in-medium modified. Due to the Dirac mass reduction with 
the baryon density, the in-medium baryon-baryon cross section is substantially reduced in 
nuclear medium. In the present work, we have extended the method of Ref. [l4] by applying 
Eqs.f l35p . fl5o]) also for the in-medium modification of the channel BiB 2 — > -B3-B4, i.e. for the 
elastic baryon-baryon collisions or resonance excitation in baryon-baryon collisions without 
outgoing mesons. 

Meson-baryon cross sections were kept always as the vacuum ones. We believe that this 
is a reasonable assumption, since the modification factor will be proportional to (m*) 2 in 
this case, while the modification factor for the baryon-baryon cross sections is proportional 
to (m*) 4 (c.f. Eq.fl3§D). Thus, the meson-baryon cross sections are less subject to the 
in-medium modifications. On the other hand, the implementation of the in-medium meson- 
baryon cross sections is more difficult, since in this case also the resonance decay widths 
should be consistently modified to preserve detailed balance in the channel MB <-> R. 

For brevity, the first scheme will be refered to as the calculation with vacuum cross 
sections below. The second scheme will be called the calculation with the in-medium cross 
sections. 



C. Three-body collisions 



Little is known about the three-body forces even in ground state nuclear matter. The 
problem of particle production in three-body collisions is even harder. There are no ex- 
perimental data on this subject and the corresponding matrix elements are not obtainable 
from data 42|. Thus, in simulating the three-body collisions, we apply a simple geometrical 



method similar to that of Refs. 
reader to Ref. [24]. 



19 



22 



23 



24J. For the details of derivation we refer the 



The geometrical method of Refs. 19 



22 



23 



24j is based on the hard-sphere collision 



picture, i.e. the potential acting between colliding nucleons 1 and 2 is assumed to be infinitely 
repulsive at relative distances d\ 2 < -Rmax and zero at d 12 > -Rmax, where -R max is equal to 
the sum of the matter radii of the colliding particles. The quantity -R max can be expressed 
in terms of the total (in-medium) interaction cross section a\ 2 of the particles 1 and 2 as 
follows: 

#max = sf^hh ■ ( 42 ) 
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The three-body collision is assumed to happen, if (i) the two nucleons 1 and 2 are about 
to collide according to the geometrical collision criterion [l| and (ii) the third particle is 
found in a sphere of radius -R max centered at the collision point of 1 and 2. The geometrical 
collision criterion for selecting the colliding pair 1 and 2 means that these particles approach 
their minimum separation distance during the given time step, and this distance is less than 
-Rmax- Thus, the three-body collision takes place when all the three particles are found 
simultaneously in the interaction volume which is here the sphere of radius -R max centered 
at the cm. of 1 and 2 in coordinate space. 

The hard-sphere collision picture reflects, in a natural way, the short-range character 
of the nuclear forces acting between hadrons. However, in a high-energy nucleus-nucleus 
collision case, one must modify this picture to account for the Lorentz contraction of the 
matter radii of the colliding hadrons. The relativistic gas of the Lorentz-contracted hard 



spheres in thermal equilibrium has been already considered in Refs. [43l . |44| in order to 
derive the Lorentz corrections to the excluded volume in the Van-der-Waalls equation of 
state. 

Therefore, in the present work we define the interaction volume as an axially symmetric 
ellipsoid contracted along the collision axis of 1 and 2 by the average 7-factor in their cm. 
frame 

712= ((7l + 72)/2)cml2 , (43) 

where 71 = pf 1 jm\ and 72 = plf/m^,. The ellipsoid has a half- axis of R max /li2 along the 
particle 1 momentum in the cm. frame of 1 and 2 and a half-axis of -R max in a transverse 
direction. 

We now select the set of all particles inside the ellipsoid which are different from 1 and 2. 
In principle, all the particles from the set plus the colliding pair form a many-body colliding 
system. However, for technical reasons, we restricted ourselves in this work to a simulation 
of three-body collisions only. Therefore, we choose only one particle from the set, namely, 
the closest particle to the cm. of 1 and 2 — as a participant of a three-body collision. We 
will denote this particle as 3 below. 

Next we simulate the actual three-body collision event of the triple 1, 2 and 3. We denote 
the initial kinetic momenta of the triple as p£, P2 and P3. Following Ref. 22j, the momenta 
of the triple are, first, redistributed microcanonically. This is done by sampling the new 
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kinetic momenta p^,, p^/ and p^/ according to the probability 

dV oc d$ 3 (pj +p* 2 +pI; p* v , p%, pi) , (44) 

where gM>3 is the three-body phase space volume element (Eq. (l38l) ). It is assumed here that 
the particles 1,2 and 3 keep their identity In particular, they stay on their initial Dirac 
mass shells: (p*,) 2 = (p*) 2 = (m*) 2 , i=l,2,3. 

After the redistribution of the momenta, the two-body collision of the particles 1 and 
2 with four-momenta p\, and p\, is simulated in a usual way. This can lead to either 
elastic or inelastic scattering, including multiple particle production through the FRITIOF 
mechanism. In Appendix A, the formal expressions for the three-body collision integral and 
the corresponding matrix element squared are given which reflects the procedure discussed 
above in the case of elastic collisions of identical fermions. 

In the actual numerical simulations, we performed the search for the third particle, which 
can be either a baryon or a meson, only if the particles 1 and 2 are both baryons or a meson 
and a baryon. For the meson-meson collisions, the search for the third particle has not been 
done. The meson-meson collisions are, however, relatively soft and can not influence much 
the high-p f part of the spectra, which are of the primary interest in our present study. 

To save the CPU time, we have also switched-off the Pauli blocking of the final state 
in all collision and resonance decay processes. We have checked by direct calculation, that 
above = 2 A GeV the Pauli blocking is negligibly small. 

III. NUMERICAL RESULTS 

We have performed the calculations for central Au+Au collisions at beam energies of 
2-20 A GeV and central Pb+Pb collisions at 30 and 40 A GeV. The time evolution of the 
systems in the cm. frame of the colliding nuclei has been followed up to 30 fm/c using a 
variable time step. The size of the time step was adjusted to reduce the spurious effect of 
multiple scatterings of the same particle within the given time step. 

In order to see an influence of the various physical ingredients of our model, four types of 
calculation have been done: (i) pure binary cascade calculation without mean field and using 
the vacuum cross sections; (ii) the calculation with the RMF, with only binary collisions 
and vacuum cross sections; (iii) the same as (ii) plus the three-body collisions; (iv) the same 
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as (iii), but with the in- medium baryon-baryon cross sections. 

It is necessary to point out, that when calculating the radius -R max of Eq. fj42|) we have 
used the cross section a\ 2 somewhat different from the actual total two-body cross section 
implemented in the model. Namely, in calculations with vacuum cross sections (i),(ii) and 
(iii), <7* 2 was set to 40 mb for a baryon-baryon collision, which is an asymptotic high energy 
value of the total pp cross section. In the case of the calculation with in-medium baryon- 
baryon cross sections (iv), we used for a* 2 the in-medium pp cross section determined by 
summation of all possible final state channels. The in-medium cross section of each final 
channel was obtained according to Eq.( l35l) . For the meson-baryon collisions, the constant 
cross section <j\ 2 = 20 mb was always used in Eq. (I42|) . This value is close to the asymptotic 
high energy n + p total cross section. The rescaling factors due to the leading quark numbers 
(c.f. Sect. II) have not been taken into account in the calculation of the interaction volume. 
To avoid misunderstanding, we stress that these simplifying assumptions have been made 
when calculating the interaction volume only, but not in determination of the two-body 
collision partners by the geometrical collision criterion. 

A. Time evolution 

Fig. [TJ shows the time dependencies of the central baryon and meson densities and of 
the total collision frequency for the Au+Au system at 20 A GeV and b=0 fm. The baryon 
and meson densities have been computed in the central lfmxlfmxlfm cube. The total 
collision frequency N tot includes the two- and (when switched-on) the three-body collisions 
and has been determined in the larger 3fmx3fmx3fm central cube to reduce the statistical 
fluctuations. The quantity N tot reaches its maximum at about 4 fm/c, when also the central 
baryon and meson densities are maximal, and drops rapidly later on. We observe, that the 
RMF reduces the maximum baryon and meson densities and also leads to a faster expansion 
of the compressed system. The three-body collisions do not influence the central densities 
and influence the total collision frequency only weakly. The weakness of the dependence of 
the total collision frequency on the three-body collisions is due to the fact, that the inclusion 
of the three-body collisions reduces the number of the pure two-body collisions, i.e. collisions 
where there is no other particles in the interaction volume of the primary colliding pair 1 
and 2. Indeed only about 10-20% of all collision events are now pure two-body (c.f. Fig. [2] 
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below) . 

It is interesting, that the inclusion of three-body collisions even reduces iVtot slightly. 
This happens, since the particle 3' and the particles emerging from the interaction of 1' and 
2' are final state particles of a three-body collision event (see Sect. IIC). The final state 
particles are not allowed to rescatter on each other before at least one of them will rescatter 
on another particle. This leads to an overall reduction of the total collision frequency, when 
the three-body collisions are included. Indeed, without the three-body collisions, the particle 
3 would be allowed to rescatter on the collision products of 1 and 2. 

Finally, using the in-medium baryon-baryon cross sections influences the baryon density 
rather weakly. However, the meson production and the total collision frequency are strongly 
decreased in this case. 

Fig. [2] presents the ratio of the frequency of three-body collisions to the total collision 
frequency (two- plus three-body) as a function of time for the Au+Au central collision at 
20 A GeV (bottom left and right panels). In order to understand this ratio better, we also 
show in Fig. [2] the time dependencies of the central baryon and meson densities along with 
their sum and the average 7-factor of the two-body collisions (top left and right panels). 
The cross section a\ 2 used to compute the radius P max of Eq. fH2|) averaged over colliding 
pairs is also plotted vs time in Fig. [2] (middle left and right panels). Calculations both with 
vacuum (left column) and in-medium (right column) cross sections are shown in Fig. [2j We 
see that the ratio N^/N to t reaches the maximum value of 0.9 for the vacuum cross sections 
and 0.8 for the in-medium cross sections. 

The ratio N^/N t ot can be estimated on the basis of the Poissonian distribution for the 
probability to find n = 0, 1,2, ... particles in the interaction volume of the colliding pair 1 
and 2 (c.f. Ref. [20]): 

P n = exp(-A), (45) 
n\ 

where A = p tot < V\i > with < Vi 2 > being the averaged interaction volume (c.f. Eq. flA2j) ): 

< V 12 >= -7T (^f^) 3/2 < 7i2 ■ (46) 

For the meson-free matter, the quantity A is proportional to the gas parameter (Eq.([T])). 
Since, by definition (see Sect. IIC), the three-body collision happens if n > 1, the following 
estimate can be done: 

N 3 /N tot ~ 1 - P = 1 - exp(-A) . (47) 
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In the nonrelativistic^ limit for the meson-free matter with A < 1, Eg. (1471) is identical to 



the result of Ref. 



19j |. The estimate (j47j) is depicted by the dashed lines in the bottom 



left and right panels of Fig. [2J We observe, that Eq. (1471) reproduces the overall behaviour 
of the directly computed ratio N 3 /N tot . In particular, one can see, that the flat maximum 
of N 3 /N tot calculated with the vacuum cross sections is caused by larger < <j\ 2 > at the 
initial stage of collision. In the case of the in-medium baryon-baryon cross sections, the 
value of < (T* 2 > drops quickly at the beginning, reaching the minimum at about 2 fm/c. 
This reduces the ratio N 3 /N tot at the initial stage of the collision and leads to the peak of 
the ratio at about 5 fm/c. At t > 10 fm/c, the meson-nucleon collisions dominate. Thus, 
< o\ 2 > is close to 20 mb in both calculations, with vacuum and in-medium baryon-baryon 
cross sections. As a consequence, at t > 10 fm/c the ratio N 3 /N tot , practically, does not 
depend on the baryon-baryon cross sections. 

Another interesting feature is that the average 7-factor (c.f. dotted lines in the top left 
and right panels of Fig. [2]) drops from initial value of 3.6 to a rather low value of 1.5 within 
5 fm/c. This reflects the transition of the fast relative motion of colliding nuclei to a slower 
thermal motion of nearly equilibrated hadronic matter. The Lorentz contraction of the 
interaction volume (Eqs. (IA2I) and PSI) ) plays only a moderate role for the equilibrated 
hadronic medium. However, at the initial nonequilibrium stage, the Lorentz contraction 
suppresses the ratio N 3 /N tot very strongly. Thus, as it should be, the very early stage of a 
relativistic heavy ion collision can be described by the binary cascade model quite well. 

Fig. [3] shows the beam energy dependence of the maximum central baryon and total 
densities and of the maximum ratio N 3 /N to t reached in central Au+Au collisions. RMF 
reduces the maximum baryon and total densities quite substantially (c.f. also Fig. []]). The 
three-body collisions do not influence these two observables. In calculation with the vacuum 
baryon-baryon cross sections, the ratio N 3 /N tot stays almost constant ~ 0.9 in the beam 
energy range from 2 to 20 A GeV. Using the in-medium baryon-baryon cross sections leads to 
a dropping N 3 /N tot towards smaller beam energies. Indeed, at smaller -Eiab, the in-medium 
reduction of the baryon-baryon cross sections is better visible, due to the smaller number of 
the meson-baryon collisions. 
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B. Comparison with experiment 



First, we address the stopping power of nuclear matter. Fig. H] shows the proton rapidity 
distributions for the central Au+Au collisions at 10.7 A GeV (upper panel) and for the 
central Pb+Pb collisions at 40 A GeV (lower panel). The cascade calculation (i) clearly 
produces too much stopping. The same efect has also been observed in earlier GiBUU 
calculations 14]. 

Including the RMF reduces the stopping power. This brings the calculation (ii) into 
closer agreement with the midrapidity proton yields. It is interesting, that, at 40 A GeV, 
the proton rapidity distribution even develops a minimum at y = indicating an onset of the 
transparency. Taking into account three-body collisions (iii) increases the stopping power 
strongly, which again results in an overestimation of the midrapidity proton yields. Ulti- 
mately, using in-medium cross sections (iv) reduces the stopping power in a good agreement 
with the data. 

The in-medium reduced baryon-baryon cross sections reduce the stopping power due to 
less thermalization, while the three-body collisions act in the opposite direction leading to 
more thermalization. Indeed, a third particle found in the vicinity of the colliding pair 
interchanges its energy and momentum with the pair: we simulated this effect by micro- 
canonical sampling of the three-body phase space (c.f. Eq. fj44l) . As a result, the relative 
momentum of the two colliding particles changes its direction according to the isotropic dis- 
tribution. Therefore, the outgoing particles are also produced isotropically in the three-body 
cm. frame. This has a strong effect on the stopping power, since the particle production 
in hadron-hadron collisions is forward-backward peaked in the cm. frame of the colliding 
particles at high invariant energies. 

The effect of three-body collisions can be even better demonstrated by plotting the trans- 
verse mass spectra of the produced mesons. Figs. [5] and [6] show the m_i_-spectra of pions, 
kaons and antikaons at midrapidity produced in central Au+Au collisions at 10.7 A GeV 
and central Pb+Pb collisions at 40 A GeV. The cascade calculation (i) overestimates the 
meson yields at midrapidity and also produces too soft K + m_|_-spectra. Including the RMF 
into propagation of the baryons (ii) reduces the yields somewhat, but does not change much 
the slopes. The three-body collisions (iii) make the m_|_-spectra considerably harder. This 
is again caused by the isotropic emission of the produced mesons in the cm. frame of the 
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colliding triple. Applying the in-medium reduced baryon-baryon cross sections (iv) reduces 
the meson yields without strong changes of the spectra shapes. 

The last fact is unexpected, since reducing the cross section should also reduce the inter- 
action volume, where the third particle is looked for, and, hence, reduce the relative fraction 
of the three-body collisions with respect to two-body collisions. However, this is true only 
for the baryon-baryon collisions, because we used the in-medium cross sections in this case 
only. Therefore, the reduction of the number of the three-body collisions of the type baryon- 
baryon- (baryon or meson) is compensated by the increase of the number of the collisions of 
the type meson-baryon- (baryon or meson). As a result, the slopes of the m^-spectra did 
not get softer after application of the in-medium baryon-baryon cross sections. 

In the case of Au+Au collisions at 10.7 A GeV, the calculation with three-body collisions 
and vacuum cross sections (iii) provides the best agreement with the experimental transverse 
mass spectra of ir + , K + and K~ . For the Pb+Pb system at 40 A GeV, the three-body 
collisions combined with the in-medium cross sections (iv) produce the best description of 
the data on m_i_-spectra of 7r~, K + and K~ . 

In Fig. [7] we present the inverse slope parameter T of the K + transverse mass spectrum 



at midrapidity vs the beam energy. To obtain T, following Refs 
midrapidity transverse mass spectrum by an exponential function: 

d 2 n 



4S|49|,|50J, we fitted the 



aexp{-m ± /T} . (48) 



m±dm±dy 

Without the three-body collisions, we underpredict the inverse slope parameter T by about 
30%. Including the three-body collisions leads to the much better agreement with experi- 
ment, except for the points at 5.93 and 7.94 A GeV, where we still underpredict the exper- 
imental inverse slope parameter by about 20%. 

Fig. [H] shows the midrapidity yields of tt + , K + , (A + S°) and K~ vs the beam energy for 
central Au+Au collisions at 1.96, 4.00, 5.93, 7.94, 10.7 A GeV and 20 A GeV and for central 



Pb+Pb collisions at 30 and 40 A GeV. The corresponding experimental data on 



antikaon and hyperon production were taken from Refs. 47J, |48|, |49|, |50j, |5l|, [52 , 



pion, 



53 



£aon, 



54 



55|. 



The calculations for the Au+Au system were done for the impact parameter range b < 3.5 



fm (5% of the geometrical cross section, c.f. 47]). For the Pb+Pb system, we have chosen 



a slightly larger impact parameter range b < 4 fm (7% of the geometrical cross section, c.f. 

H). 
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We observe that the pure cascade calculation (i) overestimates the meson and hyperon 
production. Using the RMF (ii) reduces the midrapidity meson yields by about 15%. The 
midrapidity hyperon yield is reduced stronger, by about 30%. This reflects the behaviour 
of the proton midrapidity yield shown in Fig. HJ since the mean field acts on the hyperons 
too. The introduction of the three-body collisions (iii) influences the midrapidity yields of 
the produced particles rather weakly. Finally, the in-medium cross sections (iv) reduce the 
particle production quite strongly: mesons — by about 30%, and hyperons — by about 
50%. As a result the calculation (iv) turns out to be in a good agreement with the data on 
pion and K~ production, while it underestimates the K + and hyperon yields below 40 A 
GeV. 

Fig. shows the K + /ir + ratio at midrapidity vs the beam energy. It is interesting, that 
the three-body collisions reduce the ratio quite strongly. This is due to combination of the 
two small effects visible in Fig. [HJ increase of the pion yield and decrease of the kaon yield 
by the three-body collisions. The calculation in the RMF mode with vacuum cross sections 
and three-body collisions (iii) is in the best agreement with the experimental data below 40 
A GeV. However, we fail to describe the reduction of the K + /n + ratio above 30 A GeV. 



IV. SUMMARY AND DISCUSSION 



We studied the influence of several many-body effects on particle production in heavy 
ion collisions at 2-40 A GeV. The calculations were done in the framework of the GiBUU 
model 

First, we implemented the relativistic mean field of the nonlinear Walecka model NL2 
in GiBUU. The RMF reduces the stopping power of colliding nuclei and the meson 
production. In a calculation with RMF, a part of the kinetic cm. energy of the colliding 
nuclei transformes into the build-up of the strongly repulsive time component of the tu-field. 
This leads to less compression and less entropy production by particle-particle collisions. As 
a consequence, the nuclear matter becomes more transparent when the mean field is taken 
into account (c.f. Fig.H]). This is an interesting result, since in microscopic transport models 
the mean field has been usually not taken into account (or switched-off) at high energies 



fc.f. Refs. 



m 



Second, we implemented the three-body collisions in the model by adopting the geomet- 
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rical method of Refs. 
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2J] modified to account for the relativistic effects (see 



Sect.IIC). The three-body collisions increase the inverse slope parameter of the K + trans- 
verse mass spectra to a quite good agreement (within 10%) with experimental data, except 
for the points at 5.93 and 7.94 A GeV. The additonal transverse momentum is generated 
due to isotropic emission of the produced particles in the three-body cm. frame. 

The RMF model serves also as a natural base for the construction of in-medium modified 



baryon-baryon cross sections. Here, we continued a study started in Ref. 14j. Assuming 



the matrix element of the meson production in a baryon-baryon collision to be not modified 
by the nuclear medium, we considered the in-medium modifications of phase space and flux 
factors. We also took into account the in-medium normalization of the Dirac plane wave 
bispinor. The last effect has led to a strong in-medium reduction of the cross sections. 
However, the way we extract the in-medium cross sections is very approximate, since it 
neglects in-medium modifications of the matrix element. 

We have also studied the ir + , K + , K~ and (A + E) -hyper on yields at midrapidity at 
various beam energies (Fig. [8]). Given the ambiguity in the in-medium cross sections, the 
data are reasonably well described by calculations with the three-body collisions. The same 
is valid for the K + /ir + ratio at midrapidity plotted vs the beam energy (Fig. [9]). However, 
we do not reproduce the decrease of the K + /tt + ratio above 30 A GeV. 

The problems of the microscopic transport models to describe the K + /tt + ratio have 
usually been ascribed to the excessive pion yield, while the strangeness production was well 
reproduced overall jl2,[l^, Q]- This conclusion was based on calculations within the cascade 
mode using vacuum cross sections. In the present work, by using the in-medium reduced 
cross sections, we have well reproduced the 7r + midrapidity yields, while the K + midrapidity 
yields are now underpredicted. Since pions represent the major contribution to particle 
production, we believe that the in-medium cross sections, nevertheless, provide a reasonable 
base for description of the energy transfer to the inelastic channels. The K + yield, thus, can 
be enhanced by more detailed elaboration on the strangeness production channels: BB — ► 
BYK channel below FRITIOF threshold, larger phenomenological strangeness suppression 
factor at small y/s (c.f. Eq. (J2J)). Moreover, time dependent prehadron cross sections of 
Ref. [56J could also influence the results on strangeness production. These topics deserve, 
in our opinion, future studies within the GiBUU or similar transport approaches. 

In spite of the fundamental problems with the kinetic theory at high densities, we believe 
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that our approach combining the baryon propagation in RMF with the in-medium reduced 
cross sections and three-body collisions provides a realistic description of the HIC dynamics 
in 10 A GeV domain. It could also serve to model the pre- and after-quark-gluon phase 
formation stages of a HIC at higher energies. 
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APPENDIX A: ELASTIC THREE-BODY MATRIX ELEMENT 



It is straightforward to write down the three-body elastic collision integral which corre- 
sponds to the procedure discussed in Sect. IIC: 

rentier,*! _ f 9<PP2 * * y f Q^Pti^cmU [~T * * * \ 

i coll,3bUlJ — J (2tt) 3 J (27t) 3 7 3 IV 123 ' 1 '' 2 '' 3 ''' 

x J d$ 3 {p* +pl + p* 3 ; p*,, p*,, p*,) a\r 2 } J do\, 2 ,^ v , 2 „ 

Here <j\ 2 and v\ 2 are, respectively, the total (in-medium) interaction cross section and the 
relative velocity (Eq. fl23l) ) of 1 and 2; 

V l2 = |jr(JU)V2 (A2) 

is the interaction volume, where -R max and 712 are given by Eqs. fj42l) and respectively; 
73 = P3°/ m 3'i ^3 an d ^3 are t ne three-body phase space volume (Eq. (l37|) ) and the three- 
body phase space volume element (Eq. fl38|) ). respectively; s* 23 = (p* + p 2 + P3) 2 ; c/cr*, 2 ,^ 1 „ 2 „ 
is the elastic differential scattering cross section (Eq. (l20l) ). The quantity 

cml2 £i , , A o\ 
To ^3,^ ^3 V 1 ^ 6 ) 

(2tt) 3 73 

is the density of particles in the cm. frame of the particles 1 and 2. 
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Notice, that Eg. (lAlj) takes into account the two transitions: 123 — > 1'2'3' by the micro- 
canonical sampling the kinetic momenta of 1', 2' and 3' according to Eq. fj44|) . and the elastic 
scattering 1/2' — > 1"2" afterwards. The divisions by $3 and by done in Eq. flAlD in 

order to normalize the total transition probabilities to the unity and to the branching ratio 
of the elastic channel af,y tlc * / a\ 2l respectively. Here aff^ tlc * = / da\, 2 ,^ V i 2 i,. 

Subsututing the explicit expressions for d^ 3 (p{ + p 2 + p 3 ', p*i, p?, p%) and for da\, 2 ,^ v , 2 ,i 
obtained, respectively, from Eqs. fl38|) and (120]) into Eq. (j Aip and changing the order of 
integrations, one can transform Eq. flAll) to the form of Eq. ffTTj) : 

relastic r 

i coll,3b[ilJ - J 7^3^ Cr 12 t 'l2 K 12 J (27t) 3 7 3 I V Sl23 ' m l' ' m 2' ' m 3' J 

* / I pfe / wm* ,4> w + p5 + A - * - A - *> 

In derivation of Eq. (IA4l) . we used the relations 

(y\i 2 i = 0"i»2" ^ (A5) 

1 1' 2' = I*"2" J (^6) 

valid for the elastic scattering 1'2' — > 1"2" and the detailed balance formula 



|-^1l'2'->l"2" ~~ |-^ll"2"-»l'2' • (A7) 

Finally, we have also changed the notations as 1" — > 1' and 2" — > 2'. By putting all the 7- 
factors equal to unity and also neglecting the Pauli blocking of the final states, the loss term 
of Eq. (1A4|) is reduced to the loss term of Eq.(4.17) from Ref. 24 L The comparison of the 
gain term of our Eq. (lA4j) to the gain term of Eq.(4.17) from Ref. {24 1 is not straightforward, 
since we used an assumption of the microcanonical sampling of the particles 1', 2' and 3' 
momenta, which has not been used in Ref. 241 ] . 

By comparison Eqs. ( 1A4I) and (fTTj) we get the following expression for the invariant 
matrix element of the elastic three-body collisions: 

_ _ 3! 2! g* 12 1* 2 |vr(^ max )3( 73 ) cml2 1( ^r. , * m *)?™— (AS) 

\ M Il23-1'2'3' _ ~i 7^W^Z\I7 ^3 W S 123' m l'' m 2''^3'J~^ • l A °J 



mfm2m*,m2/m3,(2(yf)' :i (27r) 4 7i2 v cr*, 2 , 
The ratio of the 7-factors in the cm. frame of 1 and 2 can be rewritten in the explicitly 
covariant form as 

(73)cmi2 2m}m*^(pJ+^) 



712 ml (m\ p\ + m\ p\) (p\ + p\) 



(A9) 
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FIG. 1: (color online) Time evolution of the central baryon density (top panel), of the central 
meson density (middle panel) and of the total collision frequency (bottom panel) for central Au+Au 
collision at 20 A GeV. Binary cascade calculation is represented by dotted lines with open rombuses. 
RMF calculation including binary collisions only and vacuum cross sections is shown by dashed 
lines with open boxes. Results with RMF including also the three-body collisions with vacuum 
cross sections are plotted by solid lines with open pentagons. RMF calculation with the three-body 
collisions and in-medium cross sections is shown by solid lines with open circles. 
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FIG. 2: (color online) Top panels: time dependence of the total (baryon plus meson) density - 
thick solid lines, meson density — dashed lines, baryon density — thin solid lines and of the 7- 
factor (Eq. (|43p ) averaged over colliding pairs — dotted lines. Middle panels: cross section used in 
calculation of the radius i? max (Eq- (|42j> ) averaged over colliding pairs. Bottom panels: the ratio of 
the three-body and the total (two- plus three-body) collision freequencies — solid histograms and 
the estimate of Eq. (|47p — dashed lines. The calculations are done with RMF including three-body 
collisions. Left column: with vacuum cross sections. Right column: with in-medium cross sections. 
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FIG. 3: (color online) Top, middle and bottom panels, respectively: the maximum central baryon 
density, the maximum central total (i.e. baryon plus meson) density and the maximum ratio of 
the three-body collision frequency to the total (two- plus three-body) collision frequency reached 
in the central Au+Au collision vs the beam energy. Various calculations are shown with the same 
notations as in Fig. [IJ 
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FIG. 4: (color online) Proton rapidity disributions for central (b < 3.5 fm) Au+Au collisions at 
10.7 A GeV (upper panel), and central (b < 4 fm) Pb+Pb collisions at 40 A GeV (lower panel). 
The experimental data for the Au+Au system are taken from Ref. [45] and correspond to 5% most 



central events. The data for Pb+Pb are from Ref. 
in Fig. [[J 



461 ] (7% centrality). Notations are the same as 
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FIG. 5: (color online) tt + , K + and K~ transverse mass spectra at midrapidity from central (b < 3.5 
fm) Au+Au collisions at 10.7 A GeV. The rapidity region is \(y — vnn) /unn\ < 0.125, where tjnn 
is the cm. rapidity in the laboratory frame. The K - spectra are multiplied by 0.1. The data are 



from Refs. 



471, 



481 ] . Notations are the same as in Fig. [TJ 
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FIG. 6: (color online) ir~ , K + and K~ transverse mass spectra at midrapidity from central (b < 4 
fm) Pb+Pb collisions at 40 A GeV. The rapidity region is \y— hnn\ < 0.2 for ir~ and \y — unn\ < 0.1 
for K 1 * 1 . The K _ spectra are multiplied by 0.1. The data are from Ref. [49(]. Notations are the 
same as in Fig. [TJ 
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FIG. 7: (color online) Inverse slope parameter of the K + transverse mass spectrum at midrapidity 
'or central collisions of Au+Au and Pb+Pb as function of the beam energy. Data from Refs. 
48, 49, 3|. Notations are the same as in Fig. (TJ 
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FIG. 8: (color online) The yield of ir + (upper left panel), K + (upper right panel), A + S° (lower 
left panel) and K~ (lower right panel) at midrapidity as function of the beam energy for central 
collisions of Au+Au at Ei„ b < 20 A GeV and Pb+Pb at E iab = 30 and 40 A GeV. The data are 



from Refs. 



47|, 



48 



49 



50, 



51 



52, 



53 



54 



551 ] . Notations are the same as in Fig. Q] 




FIG. 9: (color online) The ratio of the midrapidity yields K + /tt + for central Au+Au and Pb+Pb 
collisions. The data are from 471. l49i. l50||. Notations are the same as in Fig. [TJ 



